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Abstract 

Introduction. Essential nonlinearity of the chemical reactions of acids and bases determines the control algorithms in 
the mode of acidification or alkalization, that is, periodic dosing of a minimum volume of acid or alkali. Such regulation 
may be ineffective, specifically, it allows insufficient or excessive concentration of the controlled substance. The article 
discusses the problem of precise regulation of the hydrogen index in mini-bioreactors. It is proposed to use a digital 
model of the acidity control system to select the concentrations of topped-up solutions, determine the regulation 
methodology, and improve accuracy. The objective of the work is the assurance of required accuracy of pH regulation 
in an in vitro mini-model of the gastrointestinal tract of a static type. 

Materials and Methods. The initial block diagram of the model included accumulators and flows. It was the base for 
the main differential equations characterizing the change in volume and acidity. To correct the acidity readings of the 
resulting solution by temperature, a static model based on the polynomial approximation of experimental data using the 
least squares method was created. The structural elements of the mathematical model were investigated in the Matlab 
Simulink application package. To validate the adequacy of the mathematical model, transient characteristics were 
determined on a real system of in vitro modeling of the artificial gastrointestinal tract of poultry. 

Results. Within the framework of this work, the authors created and analyzed a nonlinear mathematical model of pH 
changes in a bioreactor taking into account external control actions. The flows of hydrochloric acid solution, alkali 
solution and drain from the reactor were presented as elements of a differential equation describing the accumulation of 
liquid in the reactor. To improve the accuracy, the solution was modified taking into account the temperature 
dependence of the hydrogen index. A dosing mathematical model based on a regulator with alkali and acid channels 
was proposed. The data obtained made it possible to generate a combined model of the pH regulation process in the 
bioreactor. The adequacy of the solution was confirmed empirically. The models of pH regulator, regulation of the 
volume of contents in the reactor and chemical reactions were shown in the form of structural diagrams. The transients 
of a mathematical model and a real control system were compared. It was established that the transient characteristics of 
the mathematical model and the real system were identical in terms of regulation time. The relative error of regulation 
of the real system was 0.35 %, and the mathematical model — 0.1 %, which corresponded to the required accuracy of 
regulation + 0.1 pH. The influence of the studied flows on the neutralization reaction was shown in the form of graphs. 
Discussion and Conclusions. The proposed mathematical model will provide selecting optimal methods and algorithms 
for regulating acidity, which will accelerate the creation of a regulator for the nonlinear process of regulating the 


hydrogen index. In the future, these developments can be integrated into a comprehensive digital model of the entire 
artificial gastrointestinal tract of poultry to optimize control algorithms (dosing, mixing, periodicity, etc.), as well as 
approximation to objects in vivo. 


Keywords: mathematical modeling, acidity, pH, in vitro modeling, in vivo modeling, control system, acidity regulation 
algorithm, digital model. 
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AHHOTalna 

Beedenue. CyjecTBenHad HesIMHeMHOCTh XMMMYeCKHX peakUMii KMCIIOT HM OCHOBaHHM oMpeyemaeT alropHTMbl 
yiipaBsIeHHA B peyKHMe NOAKMCICHHA WIM MOAMWeIaMBaHHA, TO CCTb MepHOAMYeCKON O3alHH MHHMMasIbHOrO OOBeMa 
KMCJIOTbI HIM Weroun. Takoe perymupoBaHuve MOxeT OBITh Maslo9dPeKTHBHBIM, T. €. JOMycKaeT HeOCTaTOYHy!IO WH 
W30bITOUHYIO KOHUeCHTpalluio KOHTposMpyemoro BelljecTBa. CTaTbA MocBAuleHa MpoOsemMe TOYHOTO perysIMpoBaHHA 
BOJOPOAHOTO MoKa3aTeIa B MMHH-OuopeakTopax. IIpeqnaraetcd MCHONL30BaTb WMPpOByIO MOJ[eJIb CHCTeMBI 
ylpaBsIeHHa KHCIOTHOCTbIO JIA WOAOopa KOHUeHTpalMii MONMBaeMbIX pacTBOpOB, OlIpeseyeHHA MeTOAMKH 
peryIMpoBaHHA MW MOBbILeHHA TOUHOCTH. Liens padoTEr — obecteyeHHe TpeOyemMol TOuHOCTH perymupoBaHua pH 
B in vitro MMHH-MOJeJIM 2KeJTy [OWHO-KHMIe4YHOrO TpakTa CTaTHYeCKOrO TUM. 

Mamepuanoi u memoooi. Ucxoyuas crpykTypHad CXeMa MOJeJIM BKIKOUAeT HAKONMTeIM MW MOTOKH. Ona pescTaBiaerT 
coOoi Oa3y AIA OCHOBHBIX WudpdepeHuUMaIbHbIX ypaBHeEHHH, XapakTepv3yIOWIMX H3MeHeHHe OObBeMa HM KHCIIOTHOCTH. 
Jia KOppeKTHPOBKH TOKa3aHHii KMCMOTHOCTH pe3yIbTHpyrOllero pacTBopa MO TeMIlepaType co3qaHa cTaTHuecKkaad 
MOJesIb, OCHOBaHHad Ha MOJMHOMMAIbHOM allMpOKCHMalM 9KCIePHMeHTAJIbHbIX JJaHHbIX Me€TOJOM HaMMeHBbIWIMX 
KBagpatoB. B mpuknaqHom makete MATLAB Simulink uccneqoBaHbl CTpyKTypHbIe 3JIEMeCHTbI MaTeMaTHYeCKoH 
Moyes. Ha peambHoi cuctTeme in vitro MojjeuMpoBaHuaA UCKYCCTBCHHOMO 2KeslyJOUHO-KMINedHOTO TpakTa JOMalHel 
IITHUBI Ope eeHbI MepexOMHble XapaKTePHCTHKH JIA NOATBepxK CHUA aeKBAaTHOCTH MaTeMaTHYeCKOH MOEN. 
Pezsyismamot ucciedosanua. B pamMKkax jaHHo paOoTbI aBTOpbI co3qaiu MW MpoanaM3HpoBalM HesIMHelHyt10 
MaTeMaTHYeCKy!O MOJeIb H3MeHeHHA PH B OMOpeakTope c yYeTOM BHEINUHHX yupaBIAIOWMx BosfelicTBui. [loroKu 
pacTBopa COJAHO KUHCIOTHI, pacTBOpa WesOUH U CIMBa pacTBOpa M3 peakTOpa Ipe{CTaBJIeHbI KaK 93JIEMCHTHI 
TupPepeHunalbHOro ypaBHeHHA, OMMCbIBAaIOWerO HaKOMIeHHe %KHKOCTH B peakTope. JIIa NOBbILIeHHA TOUHOCTU 
pewieHHe opadoTamm c y4eTOM TeMMepaTypHOM 3aBHCHMOCTH BOZOpogzHoro noKa3atena. Ipeqnoxena 
MaTeMaTHUYeCKad MOJeJIb NOSHPOBAHHA Ha OCHOBE perysaATopa C KaHalaMM Weoun WU KucMOTHI. lomyaeHuble JaHHbIe 
TO3BOIMIM TeHepHpoBaTb OObeAHHEHHYIO MOJeIb Mpolecca perymupoBanua pH B Onopeaxtope. AyeKBaTHOCTB 
pelleHHua NOATBepAWIM OMBbITHDIM TyTeM. B Bue CTPyKTypHbIX CXeM [MOKa3aHbI MOJeuIM: perysatopa pH, 
perysIMpoBaHua oObeMa cofepxKHMOrO B peakTOpe H XHMHYeCKHX peakuMit. CpaBHuBaloTcaA MepexOHble IMpOLlecchl 
MaTeMaTHYecKOH MOJeM UW pealibHOM CHCTeMbI yiipaBleHus. YcTaHOBJIeHO, YTO TepexoHble XapakTepHCTUKU 
MaTeMaTHYuecKOH MOeH HU peasbHOM CHCTeMbI UJCHTHYHbI TO BpeMeHH perysMpoBaHua. OTHOCHTeIbHaAd 
MOrpeLIHOCTb perysIMpOBaHHA peasbHOM cucTembI cocTaBusia 0,35 %, a MaTemaTu4eckoH Moyenu — 0,1 %, 4To 
COOTBeETCTByeT TpebyemMoli TOUHOCTH perymupoBaHua +0,1 pH. B Buge rpapukosB oKa3aHO BIMAHMe MCCIeOBAaHHBIX 


MIOTOKOB Ha peakIM1O HeliTpalu3allun. 

O@6cystcoenue u 3aknrouenua. IIpeqnaraemad MaTeMaTHYeCcKad MOJesIb MO3BONMT MOLOOpaTh ONTHMAJIbHbIC MCTO]bI 
asIFOPHTMB! peryIMPOBaHHA KUCIIOTHOCTH, YTO YCKOPUT Co3aHUe peryiaATOpa HeMHeMHOTO Mpowecca peryiMpoBaHHuAa 
BOOpozHoro HoKa3atesa. B OyayWieM 3TH HapaOOTKH MO2%KHO HHTerpHpoBaTb B KOMIMICKCHY!O IMpoBy!lo MOJeJIb 
BCerO HCKYCCTBCHHOTO Key IOUHO-KHMICYHOTO TpakTa JOMAIIHeH UTHIbI Td ONTHMU3aIIMM aNrOpHTMOB yipaBseHHa 
(03allMu, NepeMelIMBaHHA, TepHOAMUYHOCTH U T. I1.), a TAKKe AallIIpOKCHMAallHH K OOBEKTaM in Vivo. 


Koro4ueBble CsI0Ba: MaTeMaTHYecKkoe MOJCJIMNPOBaHHe, KHCJIOTHOCTH, pH, MOJCJINPOBaHHe in vitro, MOJCJIMpOBaHHe in 
vivo, CHCTeMa yUpaBJIeHHsA, asITOPpHTM peryJINpOBaHHA KHCJIOTHOCTH, mudposasa MOJICJIb. 
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Baarogapuoctu. ABTopb! BbIpaxkaloT OarofapHOcTb 3a (PHHAaHCOBY!IO NoOfWepxXKy UCcieqOBaHHit B paMKax rpaHTa 
MuuuctepctBa HaykKH HM BBICHIero oOOpa30BaHHa Poccniicxoi Wegepaunun Ne 075-15-2022-285 or 9 urona 


2022 r. «BeTepuHapHble IpoOuoTH4eckue IpenapaTbl HallpaBsIeHHOrO MOAYIHPOBaHHA 3,OPOBbA 2KMBOTHBIX)). 


Aas uuTuposanna. Joucxoi JI.1O., JlykpaHos A.J]., Dunmnosuy B. Marematuueckas MOjesIb CHCTEMBI yIIpaBsIeHHA 
pH B in vitro MoqemM 2KeTyAOUHO-KHMIe4HOrO TpakTa JoMawiHei mtuubt. Advanced Engineering Research (Rostov-on- 
Don). 2023;23(1):95—106. https://doi.org/10.23947/2687-1653-2023-23-1-95-106 


Introduction. Since the early 1980s, the scientific community studying the physiology and functional features of 
human internal organs has come to understand that in the near future, the focus of experimental techniques will shift 
from “in vivo” practices to “in vitro” systems [1]. The development of systems for hardware simulation of physiological 
processes will be of primary relevance, including for the following reasons: reduction of costs and time for conducting 
research; possibility of simultaneous conducting of the same type of experiments; and, what is not unimportant — 
increasing the ethics of research, reducing the need for experiments on living beings. By the beginning of the XXI 
century, a whole scientific direction was formed, engaged in the development of techniques and equipment for “in 
vitro” modeling and research [1]. 

Currently, in vitro models of the gastrointestinal tract are divided into two classes: static and dynamic. Static 
installations allow you to simulate the processes of digestion when using liquid nutrient media [2]. Dynamic models, 
unlike static ones, allow modeling peristaltic movements in the gastrointestinal tract, and are used when working with 
content nutrient media. 

From the point of view of the simulated volume, the gastrointestinal tract models (regardless of the previous 
classification) are divided into micro- (up to 50 ml), mini- (from 50 to 400 ml) and macro-systems (over 400 ml) [3]. 
Microfluidic systems should be mentioned separately [4], but they go far beyond the scope of this work. 

Patent analysis makes it possible to determine that the first is a patent for an invention “In vitro model of an in vivo 
digestive tract”'. This system is dynamic and precedes the market entry of the first TIM-2-type systems [5]. The 
development of this approach occurred in the models DGM [6], TIM-2 [5], HGS [7], etc. Static models, such as: DIDGI 
[8], SIMGI [9], SHIME [10], ARCOL [11], etc., moved away from complex modeling of peristalsis, and focused on 
more accurate modeling of the actual digestive processes. One of the reasons for this is the technical complexity and 
high cost of implementing systems with dynamic content compression. Another reason was the development of 
mathematical modeling methods and microcontroller control systems, which made it possible to control the medium in 
a test tube (reactor) with high accuracy and stability [4]. 

In this paper, we consider the solution to a particular problem of mathematical modeling of the pH regulation 
process for an in vitro mini-model of the gastrointestinal tract (gastrointestinal tract) of a static type [12]. The 
importance of precise control of the acidity of the chyme is determined by rather strict requirements for the conditions 


of in vitro experiments. Ensuring plausible conditions for the course of digestive processes, as well as the requirements 


for repeatability of experiments, require automatic pH control in the bioreactor with an accuracy not worse than +0.1 
units.” At the same time, the pH range can be from slightly alkaline (7.9 pH, oral cavity) up to strongly acidic (1.3 pH, 
stomach). 

The objective of this work is to ensure the required accuracy of pH adjustment in an in vitro mini-model of the 
gastrointestinal tract (gastrointestinal tract) of a static type. 

To achieve this goal, the following tasks will be solved in the work: 

— development of mathematical model of nonlinear pH change in a bioreactor taking into account external control 
actions; 


— study of a mathematical model by simulation methods; 


'US5525305A. In vitro model of an in vivo digestive tract. URL: https://patents.z0ogle.com/patent/US5525305A/en 
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“No. 075-15-2019-1880. Veterinarnye probioticheskie preparaty napraviennogo modulirovaniya zdorov'ya zhivotnykh. URL: 
https://probioticdonstu.com/itno-2021-ru/ (In Russ.) 


Information Technology, Computer Science and Management 


http://vestnik-donstu.ru 


Advanced Engineering Research (Rostov-on-Don). 2023;23(1):95—106. eISSN 2687-1653 


— refinement of the mathematical model in order to increase its accuracy by taking into account the temperature 
dependence of the hydrogen index; 

— development of a mathematical model of the dosing process based on a two-channel competitive regulator (via 
alkali and acid channels); 

— development of a unified mathematical model of the pH regulation process in a bioreactor; 

— conducting experimental studies to confirm the adequacy of the combined mathematical model and determine the 


achieved accuracy of regulation. 


Materials and Methods. To simulate the neutralization reaction, we will use the method of accumulators and 
flows [13]. The method is based on the use of fundamental or local conservation laws (accumulators) and mathematical 
description of patterns that change the amount of accumulated feature (flows). 

The block diagram of the model is described in the figure, and contains two drives: 


(1) Volume 
acid 
atts (3) Drain 
(3) Volume 
alkali 
Substance 
(4) Flow 
siciid Q> (6) Digestive 
(5) Flow pubstaliee processes 
HC1 


Fig. 1. Block diagram of the model 
The reactor as a physical object simultaneously acts as a storage of the amount of liquid (Q1) and a storage of the 
amount of hydrogen ions (Q2). This is a container in which a controlled digestion process takes place, accompanied by 
a change in pH and a change in fluid volume. 
There are three threads acting on the first drive: 
— the first stream is the volume of acid topped up; 
— the second stream is the volume of added alkali; 
— the third stream is the volume poured out of the reactor. 
The second drive is also affected by three threads: 
— the fourth flow is the supply of sodium hydroxide by a pump [2]; 
— the fifth stream is the intake of hydrochloric acid. [2]; 
— the sixth stream is a potential change in pH during digestion. However, from a methodological point of view, it is 
a disturbing effect, and is not discussed in this paper. 
Now, let us consider the mathematical implementation of the model for describing physico-chemical processes. 
The volume of liquid in the reactor is described by the following equation: 
V=S-h. (1) 
Therefore, the change in the volume of liquid in the reactor: 
dV =S-dh. (2) 
For convenience, let us call the streams: 
— the hydrochloric acid solution flow: q1=G1; 
— the flow of the alkali solution (sodium hydroxide): q2= G2; 


— the discharge flow of the solution from the reactor: q2 = G3. 
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Thus, the differential equation describing the accumulation of liquid in the reactor: 
dh_1 
To determine the dynamics of changes in the concentration of the solution in the reactor, it is required to find the 
differential of the amount of substance according to the following formula: 
v=V-C, (4) 
where V is the volume of the substance, 1; C is the molar concentration, mol/l. 


From (4) and (1), we obtain: 


# -£(s-h-c)=s(cBrnZ), (5) 
dt dt dt dt 


The change in the concentration of the solution is influenced by the flows of alkali and acid with a certain 


concentration, then: 
dh dC 
s(c-Ben E) = (6,-Cyou -Gs-Gra) (6) 


From equality (6), we express the change in concentration and simplify the expression: 


dC 1 
dt = rare (Cynon =) 4G (C=Caa)) . (7) 

Thus, the system of differential equations describing the dynamics of changes in the volume and concentration of 
the solution in the bioreactor is presented in formula 8: 


dh_ 1 


a = 5G +G, —G,), 

(8) 
dC 1 
FG Cul Canon ~C)+ G:C Cc). 


To obtain the concentration of H ((H+]) cations from the solution of the system of equations (8), let us recall the 
basics of chemistry. The concentrations of hydrogen ions ([H+]) and hydroxide ions ([OH-]) in distilled water at 25 °C 
are equal and amount to 10-7 mol/l, this directly follows from the definition of the ionic product of water, which is 
presented below [11]: 

K = [H*] x [OH] =10'mol’/?? (at 25 °C). (9) 

It is also worth considering that: 

T =[OH ]-[H"]=[Na]{Cr]. (10) 

Now we find [H+] (denoted as Cy) through the square root of the previous expressions, taking into account the 


concentration of the solution according to equation (4) of our mathematical model [14]: 


C,, =VK,T =0, (1) 


Where ions i = 4. 
To create a control system with feedback on acidity in the form of a hydrogen index (pH), it is necessary to calculate 


a negative decimal logarithm from the concentration calculated according to equations (7). Thus, we get the acidity of 


the solution under normal conditions, that is, at 25 °C. 
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But the contents of the bioreactor are not always in normal conditions. For example, to create an artificial 
environment characteristic of the in vivo environment of the gastrointestinal tract of poultry, a temperature of 42°C is 
maintained in the reactor’. Therefore, experimentally and on the basis of available research* a matrix of changes in the 
pH of buffer solutions depending on temperature was obtained on high-precision pH meters, and the readings were 
obtained using a measuring system with thermal compensation for the electrode used? [13]. 


Based on these data, a static model was constructed, it was obtained by polynomial approximation by the least 
squares method. A visual representation of the pH dependence on temperature is shown in Figure 2: [14-16]. 


Real pH 
14 


12 


10 


2 
pH 25 °C 


Fig. 2. A surface describing the dependence of pH on temperature (12) 


F(x) = p00+ pl0-x+ p0O1- y+ p20-x* + pll-x- y+ p02: y? + 12) 
+p30-x° + p21-x?- y+ pl2-x-y’. 


The coefficients of the equation are presented in Table 1. 
Table | 


Coefficients of the equation (11) 


Coefficient Value 
poo 0.154 
plo -0.00511 
pol 0.899 
p20 -2.998 - 10° 
pli 0.00333 
po2 0.0130 
p30 -3.226 - 107 
p21 1.588 - 107 
pl2 -4.98. 107 


To determine the transients in the system, we will build a mathematical model in the MATLAB Simulink 
application package (Fig. 3): 


> Automatic Temperature Compensation in pH Measurement. _ https://www.horiba.com/esp/water-quality/support/technical-tips/bench- 
meters/automatic-temperature-compensation-in-ph-measurement/?utm_source=uhw&utm_medium=301&utm_campaign=uhw-redirect 

Toa 

Ibid. 


3 


Gnaiger E. pH measurement and temperature dependence of pH. URL: 
https://www. bioblast.at/images/archive/b/be/20190329092613!MiPNet08.16_pH-Calibration.pdf 
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V_ pH(25C) 


Static model pf pH 
adjustment from 
solution temperature 


Model of the dosing 
control system 


Volume 
change model 


Neutralization 
reaction model 


Fig. 3. General structure of the mathematical model in MATLAB Simulink 


Let us consider the created system by blocks. The first block is a block of a two-channel acidity regulator with 
feedback. Noise has been introduced into the feedback circuit to simulate random events of a real system for measuring 
acidity and measuring instrument errors. 


kite 


Switch H Discrete Controller H 


Regulation error 


Real pH 


Real pH GC) 


Fig. 4. Structure of the pH regulator, block “Model of the dosing control system” 


The second unit is a system for monitoring the volume of solution in the reactor. The “Drain volume controller” 
function simulates the activation of pumping. If the volume of the solution approaches the maximum volume of the 
reactor, the function gradually resets the solution to the required volume. 


x 1/S (reactor) 


Volume to h 


Fig. 5. Structure of the reactor content volume control model, block “Volume change model” 
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The third block is the main one. It implements the calculation of the chemical reaction of acid and alkali when added 
to the reactor and calculates the pH of the final solution under normal conditions. The function called “Neutralization” 
calculates the concentration of hydrogen ions in a solution [12]. And the “Concentration in pH” block implements the 


conversion of concentration into a hydrogen pH indicator. 


pH (25 C) 


Concentration in pH 


-( Neutralization) 


Fig. 6. Structure of the model describing the neutralization reaction 


The fourth block enables to take into account the temperature changes of the solution with the function of the 


statistical model in the “Temperature correction” block. 


Temperature correction 


Fig. 7. Structure of the model describing the change in acidity from temperature, block “Static model of pH adjustment from solution 


temperature” 


To determine the adequacy of the transient characteristics of the resulting mathematical model, experiments were 
conducted on a simulation model. The real model is represented by an automated control system for the artificial 


gastrointestinal tract of poultry with one control unit and a bioreactor. Figure 8 shows the experimental equipment and 


its description. 
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a) b) 
Fig. 8. Research equipment and its description: 
a — bioreactor environment control units; b — elements of bioreactor. 1 — controlled agitators; 2 — heaters; 
3 — ports with sensors and a system for topping/draining solutions 


The methodology of the experiment is as follows. The mathematical model takes into account the characteristics of a 
real system, such as: 

— set acidity; 

— reactor parameters (base area and height of the bioreactor); 

— concentrations of added substances (acids and alkalis); 

— initial concentration of the solution; 

— initial volume of the solution; 

— the temperature of the solution throughout the experiment; 

— discreteness of the control process [14]. 

The mathematical model is an ideal system in which mixing of solutions in the reactor occurs instantly and 
uniformly. In real conditions, this is impossible, therefore, forced mixing is needed for the uniform dissolution of the 
added acid or alkali. In the simulation model, this problem is solved by a magnetic stirrer controlled via a SCADA 
system. 

Unfortunately, this type of mixing device induces electromagnetic interference on the pH electrode. To eliminate 
interference, the mixing device was switched off for a short time, and the discreteness of measurements and the 
operation of the regulator was 20 seconds, which is due to the experimentally established time of normalization of the 
pH readings of the electrode in the reactor (according to the documentation of the electrode, sampling can reach 1-2 
minutes).°’ The regulator of the real system is limited in the maximum volume of the dosed acid or alkali in order to 
ensure safety, which was taken into account in the model regulator. 

Experimental Results. Figure 10 shows a comparison of the transition process of a mathematical model and a real 
control system. The pH electrode is a sensitive element that is susceptible to receiving external electromagnetic 
interference. 

One of the experiments was carried out at a given acidity of the solution in the range of 2.15 + 0.05 pH. The acidity 
of the initial solution was 6.2 pH. The temperature of the solution in the reactor: 22.25 °C. Solutions of 0.1 mol/liter 
NaOH and 0.1 mol/liter HCl were used to change the acidity. The dosing algorithm consisted in discrete proportional 
regulation. 


° ph-4502c pH meter calibration notes. URL: https://tlfong01.blog/2019/04/26/ph-4502c-ph-meter-calibration-notes/ (accessed: 30.11.2022). 
7 EZO™ pH Circuit Datasheet. URL:https://files.atlas-scientific.com/pH_EZO_Datasheet.pdf (accessed: 30.11.2022). 
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0 100 200 300 400 500 600 Time, s 


Fig. 8. Transient characteristics of a mathematical model and a real system when regulating from 6.2 to 2.15 pH 


Figure 11 shows the flow integrals from differential equation (7). 


6-10° T T T T x T T T T , — GOH] 
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bp an --- G, [fluid] 
E 
3 
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amy 
4 4 = 


350 400 450 ‘Time, s 
Fig. 9. The effect of each flow on the neutralization reaction during experiments 
Discussion and Conclusions. Analysis of the simulation results have shown that flows with a solution 
concentration in the reactor do not make a significant contribution to the regulation of the entire system, therefore, they 
can be neglected for simplification. Then the system of differential equations of the digital model (8) will take the 


following form: 


h 1 
TE eG 26.26); 
d Ss : 
dC 1 aa 
rs = a *Cyaon ~ Gr *Cuer ) 
According to the results of the experiments, taking into account the concentration of the solution in the reactor, the 
relative error of the steady-state acidity of the solution for the mathematical model was 0.1 %, and for a real system 


with the same algorithm 0.35 %. 


Thus, taking into account the error of the electrode readings and the neglect of the mixing process of solutions in the 


mathematical model, the absolute error of the control results is within acceptable limits and does not exceed + 0.1 pH of 
the set value. This means that the developed mathematical model adequately describes the neutralization reaction curve, 
and in the future, it will allow choosing optimal control algorithms, which would be a laborious process on a real 


system. 


DYu Donskoy, et al. Mathematical Model of the pH Control System in an In Vitro Model of the Gastrointestinal Tract of Poultry 


The developed model solves the problem and is a mathematical tool for the development of methods and algorithms 
for pH regulation, taking into account the temperature of the contents of the mini-reactor, the initial concentration of the 
main solution and the concentrations of top-up solutions, which will improve existing control algorithms. 
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